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Abstract 

A study of some decomposition schemes of the molecular energy into one- and two- 
center contributions published in the above mentioned journal highlights the impor- 
tance of a 'dilemma' raised in such decompositions. Even more, it has been recently 
assigned a prominent role in the promotion energy mechanism. This critical comment 
clarifies the validity of such a 'dilemma'. 



Two decomposition schemes of the molecular energy into one- and two-center 
contributions, E = J2a + Y1iA<b ^ab have been reported in Refs. They 
only differ in the treatment of the kinetic energy terms what has been considered as a 
dilemma [1]. One of them (model 1) yields too large numbers (in absolute value) but 
presents an expected distance behavior; the other (model 2) provides numbers on the 
chemical scale but exhibits a counterintuitive distance dependence. The aim of this 
report is to study the behavior of both schemes in a wide interval of internuclear 
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distances in order to check the existence of that dilemma. The two partitioning schemes 
are P, 
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in which the superscripts (1) and (2) stand for models 1 and 2, respectively. In these 
formulas /i, z/, p, r, .... are non-orthogonal atomic orbitals which constitute the basis set. 

= T - ^ is the one-electron term of the Hamiltonian related with the atom A, 
being T the kinetic energy operator and ^ the electron-nucleus interaction; Za is 
the nuclear charge of atom A and ta the electron-nucleus distance. D is the usual 
first-order reduced density matrix and [i/r|/ip] the two-electron integrals in the [12|12] 
convention. Notation /i G A indicates that orbital p is centered on the atom A. 

We have tested the behavior of both models studying the dissociation of the 
ethane molecule (eclipsed in its singlet ground state) into two methyl groups. The 
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calculations have been performed using the GAMESS package [3j at the same level of 
approximation (restricted Hartree-Fock (RHF)) with the same basis set (6-31G**) that 
were used for this system in Ref. [1]. The results are gathered in Fig. 1, where the 
distance dependence of the two-center energy component Eqc is reported for the two 
schemes. Both curves differ in the value of their energy minimum and in its localization. 
The partitioning of the kinetic energy into one- and two-center contributions (model 
2) leads to a curve appreciably softer than when this quantity is only treated as a one- 
center feature (model 1). Model 1 predicts an unphysical value Ecc = 447.95 kcal/mol 
(0.7139 a.u.) at the minimum of the curve which is very far from the chemical scale, 
while a more reasonable value Eqq = 129.57 kcal/mol (0.2065 a.u.) is attained within 
model 2, in agreement with the experimental binding energy value 100.5 kcal/mol 
(0.1602 a.u.) |4]. These results show that the assignment of the kinetic energy only to 
one-center contributions causes an abrupt lowering of the energy Ecc which may be 
attributed, as the author claimed, to a fast decay of the kinetic energy contributions 
in comparison with those of electrostatic nature. The two curves have been fitted by 
a Morse potential [5] showing that the points in the dissociation hypersurfaces follow 
an expected behavior in both cases. Hence, we may argue that the counterintuitive 
distance dependence of the model 2, which is the conclusion reported in Refs. IH Ej, 
arises from calculations there performed in a too short interval of distances CC situated 
in the lowering part of the curve. 
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Figure 1: Dissociation of the ethane molecule into two methyl groups in eclipsed conformation 
at RHF level of approximation using the 6-31G** basis set. The geometry of each CH3 group 
has been optimized at stretching. 
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